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Abstract 

We present a high-resolution photoluminescence study of Er-doped Si02 sensitized with Si 
nanocrystals (Si NCs). Emission bands originating from recombination of excitons confined in 
Si NCs, internal transitions within the 4f-electron core of Er 3+ ions, and a band centered at 
A ~1200 nm have been identified. Their kinetics have been investigated in detail. Based on these 
measurements, we present a comprehensive model for energy transfer mechanisms responsible for 
light generation in this system. A unique picture of energy flow between the two subsystems is 
developed, yielding truly microscopic information on the sensitization effect and its limitations. In 
particular, we show that most of the Er 3+ ions available in the system are participating in the 
energy exchange. The long standing problem of apparent loss of optical activity of majority of 
Er dopants upon sensitization with Si NCs is clarified and assigned to appearance of a very effi- 
cient energy exchange mechanism between Si NCs and Er 3+ ions. Application potential of SiC^Er 
sensitized by Si NCs is discussed in view of the newly acquired microscopic insight. 

PACS numbers: 78.55.-m, 31.70.Hq, 78.67.Bf, 73.22.La 
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I. INTRODUCTION 



SiC>2 matrix doped with Er 3+ ions and Si nanocrystals (Si NCs) is intensively investigated 
as an interesting system where room temperature 1.5 /xm Er-related emission can efficiently 
be induced by non-resonant excitation; some time ago it has been noted that photo- and 
electroluminescence of Er 3+ ions in SiC>2 can be effectively sensitized with Si quantum dots 



J ASM 



5j. In photoluminescence, upon illumination, incoming photons are predominantly 
absorbed by band-to-band transitions in Si NCs. Since the indirect band structure of Si is 
preserved also in the nano crystalline form 6j, electron- hole pairs generated in this way 
are characterized by a relatively long lifetime. This enables energy transfer to Er 3+ ions 
located in vicinity of Si NCs. In that way a channel of indirect excitation of Er dispersed 
in Si02 is created and the 1.5 /ini Er-related emission appears. Its temporal characteristics 
comprises a microsecond rise time, corresponding to Si NC-to-Er energy transfer, followed 
by predominantly radiative and temperature independent decay in the millisecond range, 
characteristic for Er 3+ ions in Si02- It has been concluded that dispersion of Si nanocrystals 
in Er-doped Si02 matrix Si02:(Er, Si NCs) combines to a certain extent positive features of 
Er-doped crystalline Si with those of Er-doped Si0 2 : 

• In contrast to the situation for Si02:Er, introduction of Si NCs enables indirect ex- 
citation of Er 3+ ions. This process is non-resonant and relatively efficient, with an 
(effective) excitation cross-section of a ~ 10~ 17 — 10 -16 cm" 2 , which represents an 



increase by a factor 10 3 



4| in comparison to Si02:Er. 



• In contrast to the situation in crystalline Si:Er, emission from Er 3+ ions does not 
suffer from thermal quenching and is readily observed at room temperature, similar 
to Si0 2 :Er ^Q. 

These promising characteristics rose considerable hopes on possible applications of the 

n h n 

Si02:(Er, Si NCs) for Si photonics and specific devices have been proposed [3, [8|, [9|. Par- 
ticularly attractive is the prospect application of Si02:(Er, Si NCs) for development of a 
flash-lamp pumped optical amplifier — a much welcome replacement for the currently used 
fiber amplifier which requires resonant and high power laser pumps for its operation. In 
order to achieve that, the Si NC- induced sensitization process of Er emission in Si02 has to 
be thoroughly understood. 
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In the past, Forster (dipole-dipole) mechanism [10| has been proposed in order to explain 
the energy transfer from Si NCs to Er 3+ Jill . and different locations of Er in respect to Si 
NC have been considered Jl3,[l^]. In addition to the relatively slow (/is range) NC- mediated 
Er 3+ excitation, also a much faster (100 ns range) and usually dominating process has been 



concluded 



15] . Its physical origin has been considered recently 16| and participation of 



special luminescence centers in erbium excitation has been postulated in particular for low 
temperature annealed samples. However, in spite of a considerable progress by both mod- 
elling and experiment 4j, [l7| , many issues still lack sufficient explanation. In this field, while 
reports on optical gain have been published [3], detailed investigations revealed important 
drawbacks of the SiC>2:(Er, Si NCs) system: it has been realized that (i) only a relati vely 
small part of all the Er 3+ ions is susceptible to the indirect excitation via Si NCs [jj], 20| . 
and (ii) upon introduction of Si NCs, a considerable portion of Er dopants loses so-called 
"optical activity" and does not contribute photons regardless of the excitation mode (via Si 
NCs or direct by resonant pumping) 19[ . Therefore it has become clear that dispersion of Si 
NCs in Si02:Er is a challenging and complex physical system, which has to be understood, 
and possibly engineered at a microscopic level, before device applications can be considered. 



In our previous contribution 



2l| . sub-microsecond Er-related luminescence from the 



Si02:(Er, Si NCs) system was reported, and an Auger-facilitated energy transfer process 
between carriers within the quantized levels of the NCs and Er 3+ was proposed as the ex- 
citation and de-excitation mechanism. It was also shown that up to ~ 50% of the total Er 
content is involved in this process and contributes to the sub-/is emission. In this paper, we 
study in depth temporal details of emission bands from Si02:(Er, Si NCs). Using optical 
excitation with nanosecond pulses and time-correlated photon counting detection mode, we 
resolve true kinetics of emissions related to Si NCs and Er 3+ ions present in the investigated 
material. Based on this information, we propose a complete microscopic scenario of energy 
transfer processes in the Si02:(Er, Si NCs) system. The chosen approach allows us to iden- 
tify emissions from most of the Er dopants available in the material. Using a theoretical 
model, we discuss microscopic aspects of energy transfer processes between Si NCs and Er 3+ 
ions in relation to their (mutual) location in the Si02 matrix. 
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II. EXPERIMENTAL 



For the purpose of this study, a series of Si NC- and Er-doped Si02 2 jum layers were 
prepared by radio frequency co-sputtering on a SiC>2 substrate. The samples were character- 
ized by different concentration of Er dopants and size and concentration of Si NCs. Details 
of sample preparation procedure can be found in Ref. Q]. Following preliminary selection, 
the detailed investigation of recombination dynamics have been conducted on a particular 
sample with the most intense 1.5 /zm Er-related emission. It was characterized by atomic 
concentrations of 0.11% of Er, 41.8% of Si, and 12.6% of excess Si. After the sputtering 
procedure, the sample underwent annealing at 1100 °C which resulted in formation of Si 
NCs with the average diameter of 3.1 nm with a size dispersion of ~ 14%, and a density of 
[NC] ~ 4.1xl0 18 cm -3 . A twin sample without Er dopant and the same concentration and 
size of Si NCs was used for the absorption measurements. 

The photoluminescence (PL) experiments were performed under pulsed excitation in the 
visible provided by a tuneable optical parametric oscillator (OPO) pumped by the third 
harmonic of a Nd:YAG laser, with pulse duration of 5 ns and repetition rate of 10 Hz. In 
the UV range experiments, the excitation was provided by the third harmonic of a Nd:YAG. 
The samples were placed on a cold finger of closed-cycle cryostat and the measurements were 
taken at 10 K and room temperature (RT). PL spectra were resolved with a TRIAX 320 
spectrometer and detected with an InGaAs photomultiplier tube (PMT) with flat response 
in the visible to near-infrared range or a Ge photodiode (Edinburgh Instruments EI-A) 
connected to a digital oscilloscope where signal integration was done. For time-resolved 
measurements of PL dynamics, the PMT was working in time-correlated photon counting 
mode with temporal resolution up to 2 ns. 

The absorption measurements were performed in two experimental configurations: Direct 
transmission measurements under OPO pulsed illumination, and in a UV absorption photo- 
spectrometer. 
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III. RESULTS 



A. Photoluminescence spectra 

Figure [T] shows the PL spectra of the investigated sample at RT and at 10 K. The 
measurements have been performed under pulsed excitation at A = 450 nm (2.75 eV), 
i.e. not in resonance with any internal transition of Er. The total time-integrated PL 
response of the sample is given, which reflects the number of photons emitted. Both, Er- 
and Si NC- related PL bands are observed simultaneously, with their mutual intensity ratio 
determined by the Er concentration, in agreement with earlier studies j^j]. In addition, at 
low temperature, a weaker band centered at A ~ 1200 nm can be seen; it has a short decay 
time constant and therefore its time-integrated intensity is small compared to the other two 
bands. 

We note that the Er-related emission band broadens at RT, keeping its wavelength in- 
tegrated intensity practically constant, also in concordance with previous reports. The 
broadening at RT is produced mainly in the higher energy side of the spectrum, and the 
full width at half maximum (FWHM) increases from ^6.7 meV to ~ 20 meV at 10 K and 
RT, respectively. This could be due to population of upper states of the multiplet 4 Ii3/2 by 
thermalization. 

The Si NCs show a broad excitonic-related PL band. If we take into account the NCs 
size (radius R = 1.55 nm) and dispersion (« 14%), the position and width of the excitonic- 
related PL band can be estimated on the basis of the calculated band gap of the NCs as a 
function of their size j^j. The emission energy is given by E = E^ 1 + E le + E^ - E, 



excit 



hu ~ 1.5 eV (with Eg 1 corresponding to the energy band gap of crystalline silicon, Ei e , 
and E excit the electron and hole leve ls q uantization corrections and the excitonic correction, 
respectively — taken from reference (22 1 — , and the phonon energy fko = 60 meV), i.e. a 
band centered at 850 nm, ranging from pa 820 nm to ~ 870 nm taking into account their 
size dispersion. Comparison with experiment shows that the center of the band is in very 
good agreement with the expected value, but the experimentally recorded band is broader 
than expected from the calculations, with the higher energy side of the band suffering from 
strong temperature quenching. Since optical excitation was provided by short wavelength 
photons of 450 nm (2.75 eV), thus creating hot carriers in the higher energy states of the Si 
NCs, then it is reasonable to consider the origin of the higher energy side of the NC-related 
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PL band as arising from "hot" carriers in the upper electron (and/or hole) levels of the Si 
NCs. 

Superimposed on the Si NCs PL band recorded at 10 K, peaks with positions at 1.26 
eV and 1.57 eV are observed. These correspond to emissions from the second and the 
third excited states of Er and disappear at room temperature. We note that in view of the 
exclusively non-resonant excitation mode, the observation of these peaks indicates that the 
Er 3+ excitation proceeds, at least partly, via higher excited states. 

The effective lifetime of the excited state of Er in SiC>2 with Si NCs lies within 2-3 ms. 
Decay of excitons in Si NCs is governed by a stretched exponential function, with a long 
final tail, with decay time constant of 20 fis to 50 fis as will be shown later. The decay 
of the band centered at A « 1200 nm is clearly faster, being predominantly non-radiative; 
the lifetime shortens further at a higher temperature, leading to a strong reduction of the 
time-integrated PL intensity of this band at RT. 



B. High-resolution photoluminescence kinetics 

Using the time-correlated photon-counting technique to resolve the time evolution of the 
PL transient revealed new insight into excitation and de-excitation dynamics in the Si02:(Er, 



Si NCs) system [2l[. This experimental technique allows for simultaneous recording of very 
different dynamical ranges without suffering from signal distortion. Panel a in Fig.[2]presents 
a time-resolved PL spectrum of the excitonic PL band from Si NCs. The contour plot shows 
how the peak of maximum intensity drifts to longer wavelengths — i.e. lower energy — at 
longer delays. Consequently, the PL decay constant of the luminescence increases for longer 
wavelengths. This is understandable when we consider the NCs size distribution: excitonic 
PL dynamics in smaller NCs is faster than in larger quantum dots. In the Si NCs under 
consideration the phonon-assisted radiative transition dominates over the direct radiative 
transition. The probability of such a transition increases rapidly with the decrease of the 



NCs size (see Ref. |22|). In the panel b of the figure, the time- resolved PL spectrum of Si 
NCs is shown in the sub-microseconds time range; we can observe the center of the band at 
A = 775 nm (1.6 eV), in contrast to the previous figure. If we compare the spectra measured 
at 100 ns and at 100 /is after the excitation pulse, a difference of ~ 200 meV separates 
the centers of the two bands. We can assign the sub-microsecond PL to recombination of 
carriers from the higher Si NC states. 
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Figure [3] shows measurements related to the broad PL band appearing in the spectral 
range between the NC- and Er-related emissions. In the inset, the spectrum recorded at 200 
ns after the excitation pulse and at T = 10 K, shows the broad band centered at A ~ 1200 
nm. The lifetime of this PL band, as mentioned before, shortens at higher temperatures, 
becoming practically negligible at RT. This can be seen in the main panel, where the RT 
PL transient kinetics recorded at the maximum of the band (1200 nm) is compared to the 
PL at 1535 nm. One can see that after 100 ns the decay of the PL at 1535 nm can be 
attributed to the Er-related PL. The origin of the band centered at A ~ 1200 nm is usually 
identified as recombination at defects, but its broadness and fast kinetics seem to contradict 
this identification; future investigations — out of the scope of this paper — will elucidate this 
point. 

In Fig. |H the high-resolution kinetics of the excitonic luminescence from Si NCs can be 
followed in the panel a. An intense PL signal with a fast decay, also characterized in the 
sub-/xs time domain, is observed after the excitation pulse. This is stretched until the final 
slower tale — of the order of 30 to 50 fis — is achieved, as shown in Fig. [21 Panel b of the 
same figure shows the Er-related PL kinetics for the first microsecond after the excitation 
pulse, at T = 10 K and RT. 

In our previous contribution 2l| we have analyzed detailed kinetics of Er-related PL 
sensitized by Si-NC and have shown that it exhibits 3 separate regimes: 
Regime I - for t ^ 1 [is: a strong emission appears immediately following the laser pulse 
and then rapidly decays towards a temporary minimum. 

Regime II - for 1 /is ^ t ^ 10 /is: Er-related PL intensity rises slowly again and reaches a 
broad maximum whose amplitude is at least an order of magnitude smaller that the initial 
value in Regime I. 

Regime III - for t > 10 /is: PL intensity slowly decays with a millisecond time constant. We 
have also shown that the microsecond rise of the Er-related PL in Regime II (r^ e ~ 3.5 ^s) 
is paralleled by a decay of Si NC PL, which can be characterized by a stretched exponent 
with tnc ~ 1-2 fis, illustrating in this way the relatively slow energy transfer between Si- 
NCs and Er 3+ ions. This decay of NC-related PL slows down once Er emission attains a 
maximum. At this stage the NC-mediated excitation of Er is completed or saturated 
and therefore further decay of the NC-related PL is decoupled from Er. Careful inspection 
revealed that the local minimum of Er PL intensity from which the "slow" rise begins, is 
equal to approximately half of the maximum intensity attained in Regime II. We recall that 
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such an initial amplitude of Er PL attained very shortly after excitation - Regime I - has been 



2l| — , we investigate the 



reported also before 23j. Based on the total evaluation of optically active Er fraction, we 
conclude that the "residual concentration" of Er 3+ ions involved in the slow excitation/slow 
recombination process - Regimes II and III - amounts to ~ 0.2% of the total Er content. 

To validate the idea of hot carriers in the Si NCs being responsible of the sub-microsecond 
Er-related luminescence — postulated in our previous contribution 
intensity ratio between the fast and the slow components of Er PL for different excitation 
wavelengths and fluxes. Both, fast (sub-microsecond) and slow (milliseconds) components 
of the Er-related transient PL must be governed by different excitation and de-excitation 
mechanisms. The intensity ratio of these components is thus likely to depend on the exci- 
tation conditions. Figure shows full, high- resolution kinetics of Er 3+ PL, under different 
excitation wavelengths and powers, normalized for the maximum intensity of the slow com- 
ponent. One can observe a higher fast-to-slow intensity ratio when higher power or larger 
photon energy quantum are used for excitation. 

C. Excitation cross section measurements 

The excitation cross section of Si NCs and Er PL, o~nc an d cr Er , respectively, is of crucial 
importance to understand the excitation processes and energy transfer between Si NCs 
and Er 3+ ions. In order to gain insight into this aspect, flux dependencies of Si NCs- and 
Er-related PL signals were recorded; a reference sample without Er doping was used for 
the NC-related PL measurements. In Fig. [6], PL flux dependencies for different excitation 
wavelengths are shown, recorded at 1535 nm and 912 nm, for the Er-doped and Er-free 
samples, respectively. The curves have been fitted to the excitation dependence of PL 
intensity, derived from the rate equations for pulsed excitation, 



L PL 



oc N* = N(l - e -^ A ') (1) 



where Ipl is the time-integrated PL intensity, in arbitrary units; N* is the number (concen- 
tration) of excited emitters (Er 3+ ions or Si NCs); N is the total number (concentration) 
of excitable Er 3+ ions or NCs; a is the effective excitation cross section; <fi is the photon 
flux of excitation; and At the laser pulse duration. Remark that formula ([1]) is valid only 



in a limit of one radiative exciton per Si NC 



191 ] . The effective excitation cross section a 



is determined from the experimentally measured dependencies depicted in Fig. |6] by fitting 
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them using formula (JTj). The resulting values, for both Si NCs and Er luminescence (ctnc 
and (Jet)i are shown in Fig. [7|as a function of the excitation wavelength. 

Directly related to the effective PL excitation cross section <Jpl, we studied the optical 
absorption coefficient a in our SiC^Si NCs reference sample. The effective (Er- and NCs- 
related) PL excitation cross section values (Fig. [7]) are plotted in Fig. [S] as a function of 
the measured Si NCs' absorption coefficient a, for each given excitation wavelength X exc . 
A linear relation between absorption and excitation can be observed in all the investigated 
range for Si NCs; in the case of Er-related luminescence, a change of linear relation is 
observed at an energy threshold, at ~ 2.6 - 2.7 eV (~ 460 - 470 nm), energy above which 
the excitation cross section grows faster than the absorption. 
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IV. THEORY 



In order to rationalize the experimental results gathered in this study, we consider the 
possible processes of energy exchange between carriers confined in Si NCs and Er 3+ ions in 
Si0 2 (outside NC). 

Energy transfer between electrons and holes in NCs and /-electrons of Er 3+ ion is imple- 
mented as an Auger process (i.e. via the Coulomb interaction). The energy can be transferred 
to an erbium ion either when a confined electron-hole pair recombines, or when an intra-band 
transition of confined carrier occurs. Both processes can be accompanied by multiphonon 
transitions to fulfill the energy conservation law, as the energy spectra of both electron sys- 
tems are discrete. The momentum conservation law plays an important role in the Auger 
recombination processes as the large momentum should be transferred by the electron-hole 
recombination due to the indirect band structure of silicon. 

Confined electrons and holes energy levels as well as their wave functions calculated 



in multiband effective mass approximation 22| are used in this consideration. Luttinger 



Hamiltonian in the spherical approximation has been used for holes and the strong anisotropy 
of the electron effective mass in silicon has been taken into account. The wave function and 



flux continuities were used as the boundary conditions [24], |25|. When calculating wave 
functions outside the NC, isotropic effective masses being equal to m and 5m were used 
for electron and hole states correspondingly. Spin-orbit splitting was neglected in both Si 
and Si02- 

The conduction band of Si has six equivalent minima in the first Brillouin zone, situated 
in the neighborhoods of the six X-points. The wave function of electron can be presented 
in the form: 

ipf, = C(r)u cl/ exp(ik „r), (2) 

where u cu and k 0u (k 0u = 0.85 x 27r/ai at ) are the Bloch amplitude and wave vector cor- 
responding to the bottom of valley v (ai at = 0.54 nm is the lattice constant of silicon). 
Envelope functions £ e (r) are found as a result of a numerical solution to the Schrodinger 
equation after separating the angular part exp(im0) (m can be any integer number) as there 
is a strong anisotropy of the electron effective mass: m\\ = 0.916m , m±_ = 0.19m . The 
electron states are sixfold degenerate for m=0 and 12 times degenerate for \m\ > 0, as two 
opposite values m = ±|m| correspond to the same energy. (This degeneracy is given with- 
out taking in account an additional spin degeneracy). So the states are marked with iVe| m | 
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where the letter e shows that it is an electron state, N is the main quantum number which 
shows the order of the energy level for given \m\. For example, the ground state is marked 
as leo, which means that this is the first electron state with m = 0. 

The approximation used in 22] leads to two types of hole states in bulk Si which corre- 



spond to twofold degenerate heavy hole band having mass m^ = 0.44m and nondegenerate 
light hole band having mass 2m/ l mz/(3m/ l — mi) — 0.1 2m (mi = 0.1 6mo for Si) 26|. The 
quantum confinement gives rise to mixing of the states. There are three types of hole states 
in spherical quantum dots: i) mixed states (hm) formed by the combination of heavy and 
light ones, ii) heavy hole (hh) states, and iii) light hole (hi) states. Each state is also char- 
acterized by the full angular momentum F (0 for the light hole states and positive integer 
for the other ones) and it is 2F + 1 times degenerate as the projection M of momentum F 
onto the quantization axis (arbitrary selected) can be any integer number having absolute 
value not larger than F. The space quantization forms a series of each type of the states 
with fixed F. So they are marked with the letters showing the type of the states with index 
indicating the value of F, and the number in a series in front of it all. For example, the hole 
state with the lowest energy is of the mixed type — lhm\. 

The calculated lower energy levels of electrons and holes confined in Si NC of diameter 
in the range 2.8 nm - 3.3 nm are shown in Fig. [91 The energy range is limited to the one of 
optical pumping used in experiments (2.85 eV). Due to the existence of relatively large energy 
separations between some neighboring space quantized levels, one phonon emission processes 
will be suppressed. This effect might lead to the so-called "phonon bottleneck" slowing 
down of carrier relaxation. And although in different nanocrystals alternative rapid cooling 
mechanism have been observed - e.g. CdSe 27| - one can expect that energy relaxation of 
"hot" carriers will in our case become slower in comparison to bulk silicon. Thus, an Auger 
excitation of erbium ions in silicon dioxide is possible, similar to the impact ionization by 
hot carriers in bulk silicon, where it plays a significant role in electroluminescence, but just 
negligibly affects the excitation of the erbium photoluminescence due to the fast energy 
relaxation of hot carriers in the bulk material. 



A. Excitation due to intra-band transition 

The probability of Auger excitation of an erbium ion situated in SiC>2 at the distance a 
from the center of a Si NC as the result of the transition of a "hot" confined carrier from 
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the state % into the state i' is given by the Fermi golden rule: 
2vr 1 



h Nf 



£ K/'> z 'l e$ l/' l )\ 2 MN^E, - E v - A ff - Nfrw ph ), 
ff 



(3) 



where $ is the potential created by the /-electron of the Er 3+ ion; /, /' enumerate the states 
of /-electrons of the ion; Aff/ is the transition energy; and Nf is the degeneracy degree of the 
/ state. The integration in the matrix elements of Eq. ([3]) is to be produced over the carriers 
confinement space and the f-electron coordinate. Due to the energy conservation law, the 
confined carrier transition is accompanied by the emission of N phonons with energy ftw^. 

In the Huang-Rhys model (the model of t wo displaced oscillators with the same fre 
quency), the phonon factor Jt{N) is given by 

1 



J T (N) = exp 



-25 N T 



exp 



N hu ph 
2 kT 



'A' 



2Sy/N T {N T +l) 



where S is the Huang-Rhys factor which in the one mode approximation is given by 



,5' 



£opt — £th 

hw ph 



(4) 



(5) 



with e opt and e t h corresponding to the optical and thermal ionization energy, respectively; 
Nt is the Bose-Einstein factor: 

1 



N 7 



exp \ —± 



(6) 



and In(x) is the modified Bessel function of order N. 

As the energy levels are highly degenerate, one should average over all the initial states 
with the energy E{ and sum over all the final states corresponding to the energy E[ in 
Eq. ([3]). All these states are actually split due to nonsphericity of NCs and other factors. 
This fact is taken into account by assuming the broadening of levels and adding the value 
SE to the argument of 5-function in Eq. (j3J) and averaging over this value in the energy 
range A# = 60 meV (we have assumed that it is about the energy of an optical phonon in 
bulk Si). In result Eq. transforms into 

f^£^£l</V,M<|e*|/ ;i ,M)|^ TW , 

a 1 M,M> ■' //' 



(7) 



where iVj is the degeneracy degree of the initial state i. M and M' enumerate the degenerate 
states of levels % and i', and final and initial energies are related through 



ph- 



is) 
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Following the Appendix, one finally gets for the probability of excitation: 



7T 



1 1 / he V e 2 



W Vi = ^=— -&[ — ) ^IiH(a)J T (N), (9) 



2v^r rad i? 4 \A ffl J e 2 2 A E ' 
where e 2 is the high-frequency dielectric constant in SiC>2 (as the energy transmitted pro- 
cesses is much larger than the phonon energies 2f|), r rad the radiative lifetime of the erbium 
ion in the first excited state ( 4 /i3/2), R is the radius of the NC, and i»'i(a) are dimension- 
less factors defined in the Appendix. From Eq. (jUJ) one gets using values r rad = 2 ms and 
R = 1.55 nm: 

W Vi = 8.3 x 10 9 f^P) Ii>i(a)J T (N) s" 1 . (10) 

It has been shown that /^(a) are actually functions of the relation a/R for holes. In the 
case of electrons, there is only a weak additional dependence on the NC size, which one can 
neglect at least in the NC size range under consideration 29J. The results of calculations of 
the factors /^(a) for electrons and holes are presented in Fig. [TU1 

The energy required to get the Er 3+ ion into the first excited state 4 Iis/2 (0.81 meV) with 
reasonable excess or shortage of energy can be covered by emission or absorption of phonons. 
The parameters of the multiphonon transition accompanying the Auger processes are not 
well defined. There is no data on existence of the electron-phonon interaction for Er 3+ ions 
in the state 4 Ii3/2 in Si02- Thus, the interaction of confined carriers with optical phonons 
should be considered. The dispersion of optical phonons in bulk silicon can be neglected 
and the multimode model of the phonon transition becomes equivalent to the one-mode 



Huang- Rhys model 26j. Phonon factor Jt{N) calculated for phonon energy equal to 60 
meV (which is about the optical phonon energy in Si) is presented in Table [J calculated 
with a reasonable value of the Huang- Rhys parameter S = 0.1 at room temperature. The 
exact value of the phonon energy and Huang- Rhys factor S are not known for the material 



considered here. We have used value S= 0.1, w 
values obtained for Er ions fluorozirconate glasses 



rich is in accordance with experimental 



Rhys parameter is shortly discussed in Ref. 31] . The interaction with optical phonons is 



30] . The influence of the value of Huang- 



forbidden for electrons in silicon. One can suppose that the interaction of confined carriers 
with oxygen vibration could be also responsible for multiphonon assisted Auger processes 
in the system under consideration. The values of phonon factor Jt(N) at hujp^ = 140 meV 
corresponding to the oxygen vibrations are shown in Table [J as well. 

In Fig.[9]the most effective excitation processes for NCs with diameter 3.1 nm are demon- 
strated. The calculated values of the probabilities for these processes as a function of the 
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distance between a Er 3+ ion and the nearest NC are presented in Fig. [TTJ 
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TABLE I: Phonon factor Jt(N) calculated with S = 0.1 for two different phonon energies fkjp^ 
and temperature T = 300 K. 

B. De-excitation of erbium by carriers confined in NCs 

When considering erbium de-excitation due to intra-band transitions of confined carriers, 
one can use the same formulae from Sec. IIV Al (and Appendix), but interchanging initial and 
final states, as well as adjusting the phonon factor Jt{N). If erbium excitation process takes 
place with carrier transition energy larger than A//' = 0.81 eV, and requires N phonons to 
be emitted, then the reverse process will be described by Jt{—N) factor, which is much less 
than Jt(N). That is why most of the processes appear either in excitation or in de-excitation 
section. 

The difference is not in the factor only. One should also take into account that the tran- 
sitions under consideration take place between degenerate states. So the total probability is 
achieved by summation over the final states and averaging over the initial ones. As the de- 
generacy degree can be different for initial and final states, the probability of excitation and 
de-excitation processes can be different even for those of them which do not need phonons. 
This difference might be considerable for confined carriers, especially for holes. Since most 
of upper levels are described by larger values of moment F, having greater degeneracy de- 
gree, the probability of erbium de-excitation should be higher for the processes which do 
not require any phonons to be emitted or absorbed. 

The most effective de-excitation processes for NCs with diameter 3.1 nm are demonstrated 
in Fig. [9] and calculated values are shown in Fig. [TU 
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C. Erbium excitation by the recombination of confined carriers 



Let us consider the excitation of Er 3+ ions by recombination of confined electron and 
hole. For the NCs under consideration (d ~ 3.1 nm) the recombination energy is larger 
than 1.5 eV. Therefore the energy transfer to the Er 3+ ion by an Auger recombination of 
such an exciton can be effective only if it causes the direct transition of the ion into the 
third excited state 4 Ig/2 (energy of transition from the ground state 4 Iis/2 is A 3 = 1.55 eV), 
the fourth 4 F 9 / 2 (transition energy A 4 = 1.9 eV) or higher excited states. One should 
notice that the Er 3+ ion can practically not be excited directly into the state 4 Ii3/2 via 
such a process, as it should be accompanied by multi-phonon emission with the number 
of phonons N ^ 5. The phonon factor Jt{N) is in this case very small (see Table [I]). In 
order to calculate the transition probability we can use formula ([3]), just assuming that 
initial and final states of confined carriers now belong to different bands and dividing the 
probability by the degeneracy of the final state since only one final state is empty if there 
is one electron- hole pair in the NC. One should also choose appropriate parameters of the 
phonon system. 

Crucial for the matrix element evaluation is keeping in mind that the value hAk of mo- 
mentum transmitted during recombination process is large. The minima of the conduction 
band in /c-space are shifted from the T point by the wave vector /c = 0.85A;x (kx is the 
Brillouin zone edge). And it was shown in Ref. 32J that the momentum to transmit is 
even larger than Hko'. hAk = 1.15hkx- Such a great momentum can only be transferred to 
the /-shell of the erbium ion by the Coulomb interaction at a distance less than the lattice 
constant of silicon. So the interaction has a contact character and is determined by the 
electron and hole wave function values at the position a of the erbium. 

Once carriers are strongly confined in the NC, and the tunnelling is weak, the interaction 
is possible either inside the NC or in its vicinity. When dealing with the Coulomb interaction 
at distances smaller than the lattice constant, no screening should be taken into account 



any more, and the effective dielectric constant value can be assumed to be e c s = 1 [26|, l32 |. 
In this case the absolute value square of the matrix element in Eq. (jBJ averaged over the 
degenerate electron and hole states can be calculated in analogy to bulk Auger processes 



33 



3J] as 



| ] ] /, z > | 2 = j </ ] ^ 2 1 /') T I <^o kc> | 2 | (a) | 2 ^ ^ k^'o 2 , (11) 



eft i M i 
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where £ e, *(r) is the electron envelope function in the initial state and for shortness of nota- 
tions the total hole wave function in the final state is written as ifjp l M (r) = ^2 m ^Mm( r ) u m 
with u m (m = —1, 0, +1) being the hole Bloch functions; = 2F' + 1 is the degeneracy of 
the hole state; |(wol u cz)| ~ 0.25 is the overlap integral between the bottom of the valence 
band F 25 , and the second conduction band with k at the position in the first Brillouin 
zone where the first conduction band has its minimum {35}]. We remark, that the quadrupole 
term plays here the main role in Coulomb interaction. Using expression f[TT|) for the right 
hand side of Eq. divided by we get the expression for the transfer probability for a 
given position of Er 3+ and radius of the NC 

W tr (a;R) = y^(— r) Q(*,R)\(uo\u cz )\ 2 ^fj T (N). (12) 

Here the factor Jjrj comes from the summation over / and averaging over /' of the absolute 
value square of the matrix element (f\z 2 \f), where r/ « 0.43 A is the radius of the 4/-shell 
of the Er 3+ ion and the unknown factor 7/ is of the order of 1. We have introduced a 
dimensionless factor Q(a; R) defined by 

Q( a ;R)= (f^ 3 ) 2 Ir^r^El&tof. (13) 

If we assume the homogeneous probability distribution for the Er 3+ ion inside the NC then 
the probability of excitation transfer averaged over the position of the Er 3+ ion inside the 
NC is given by Eq. ()12p where in place of Q(sl] R) we have 

Q™( R ) = l^m[ d 3 aQ(a;fi). (14) 

-3--^ Ja<R 

In order to calculate the average transfer probability at some distance D from the NC 
boundary we introduce 

Q S uAR) = ^JanQ(R;R), (15) 

where the integral is taken over the full solid angle Q. Then the above-mentioned probability 
is given by Eq. (fl2]) where in place of Q(a; R) we have 

Q surf ( J R) exp[-2(K e + K h )D/R]. (16) 

Here the dimensionless factors 



H 



e(h) 



'2m e{h) {U e{h) -E e{h) }Ri 



< 17 > 
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determine the decay of the electron and hole wave functions outside the NC, E e and Eh are 
the electron and hole quantization energies, m° and m° h are the electron and hole masses 
outside the NC and U e = 3.2 eV and Uh = 4.3 eV are the corresponding energy barriers at 
the NC boundary. For the considered NCs 2{k e + /%) is on the order of 10 2 (see Table ILT|) . 
Therefore the probability of the excitation transfer by the electron-hole recombination decays 
rapidly with increase of the distance between the erbium and NC. It becomes negligible at 
the distance of only several angstroms. 

For local vibrations of erbium ions in fluorozirconate glass by optical transitions from the 
higher excited states into the ground state the values around ftw^ = 60 meV and S = 0.1 
were reported in Ref. |3(|. We note that bulk optical phonons in Si also have approximately 
the same energy. These values we have used for the calculation of the phonon factor Jt{N) 
in the transfer probability (see Tabled]). For the estimation of the transfer probability we 
have also used 7^ = 1. Then Eq. (TT2]) can be written as 

W tv = 0. 8 x 10 n QJ T (N) s _1 . (18) 

We have analyzed transitions induced by electron and holes being in one of the two lowest 
states. The numerical factors Q m and Qsurf are given in Table HT1 for R = 1.55 nm together 
with the corresponding energies which should be compensated by phonons and the decay 
factor 2(k e + kh) ■ 



Transition 


Ei'i — A03 


E^i — A04 


Qin 


Qsurf 


2(k e + k h ) 


leo — > \hm\ 


-41 meV 


-391 meV 


1.34 


0.043 


98.9 


2eo — > lhmi 


186 meV 


-164 meV 


1.24 


0.098 


98.0 


leo — ► lhhx 


128 meV 


-222 meV 


1.06 


0.19 


97.7 


2eo — ► lhhi 


327 meV 


-23 meV 


0.77 


0.20 


96.8 



TABLE II: Calculated parameters of several interband transitions for NC diameter of 3. 1 nm. 

From the data presented in Table ILT1 and using Table [H Eq. (TlBl) and Eq. ( Tl8l) . one can 
see that the electron-hole recombination can effectively excite an Er 3+ ion situated inside 
the NC or at a very short distance from its boundary on a nanosecond or even shorter time 
scale. However, this ultrafast excitation process does not lead to an immediate excitation 
of the first 4 Ii3/2 excited state of the Er 3+ ion relevant for the 1.55 /im emission. Transition 
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to this state can occur only via a subsequent multiphonon relaxation process which takes 
place on a microsecond time scale. 



D. Dipole-dipole contribution for the inter-band transitions 

When an erbium ion is situated at a distance a ^> R one can expand the Coulomb 
interaction between the recombining electron-hole pair in the NC and the erbium ion in 
Si02 into the series over the coordinate of the confined carrier. In the leading order such 
an approximation results into the dipole-dipole interaction being the reason for the so- 
called Forster mechanism of excitation [u]] . In this case the large excess momentum of the 
confined pair is transferred to the boundary of the NC in the process of recombination, 
or the recombination is accompanied by emission of a phonon as in the radiative exciton 
recombination 22]. 



The probability of the excitation governed by the dipole-dipole interaction can be pre- 
sented as 

w « = Tm£-£&Z%MN,), (19) 

where Jt{Nj) is the phonon factor and d ex is the dipole momentum of the confined exciton, 
which can be estimated by using its relation to the confined exciton radiative lifetime 

1 = 4e 2 E e 3 xC^eff , x 

r- d 3fr 4 c 3 ' 1 ' 

Here _E GX is the exciton energy and the effective refraction index n e s is determined by the 



formula from Ref. 



36] 



^ eff = ( — ) ej'\ (21) 



where e^i and e m are the dielectric constants of silicon and medium, respectively, and e c g = 
(esi + 2e m )/3. 

Matrix elements doj in ffl9|) correspond to the transitions in the /-shell of Er +3 , and can 
be expressed via the corresponding oscillator strengths P^f 

d% = 9 3 f Poj, (22) 
J 2m A 0j n m 



where, in the simplest approximation, n m is the refraction index of the medium [37|. To our 
knowledge, there is no data in the literature concerning the oscillator strengths of transitions 
between the levels of the Er 3+ ion in the considered inhomogeneous media. However, we can 
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estimate them via the oscillator strengths found for several glasses and solutions 37|, |38| : 



-P03 = 1 — 3 x 10~ 7 for transition 4 I 9 / 2 — > 4 Iis/2, -P04 ~2x for transition 4 F 9 / 2 ^ 4 Iis/2- 
Based on these data we estimate 4 = 1 x 10~ 22 cm 2 , dj 4 = 7 x 1CT 22 cm 2 . 
The calculation leads to the estimation 

"fc* 1 "- 1 (£)';=■ < 23 > 

Thus, one can see that the Forster mechanism does not work effectively for the considered 
system, especially at some distance from the NC, because the radiative recombination of 
confined carriers is a faster process. Again the excitation of the first excited state of the 
Er 3+ is additionally delayed by the multiphonon relaxation from the higher excited states. 



V. DISCUSSION 

In the previous section we have shown that the Si-NC mediated excitation of Er 3+ ions 
can proceed by a variety of physical processes. For all of these, efficiency depends on the 
distance between Si-NC and Er. Therefore we have generated a simulation of the Er 3+ 
distribution as a function of the distance from the nearest Si NC — see Fig. [121 based on 
sample characteristics (Er and NC concentration, NC size), assuming a random distribution 
of both, Si NCs and Er 3+ ions in the Si02 matrix. As can be seen, the vast majority of 
the Er 3+ content is positioned outside the NCs, with only 6.6% of the Er 3+ contained inside 
the NCs volume. In addition to this statistical prediction, we note that during crystalliza- 
tion process a considerable part of the Er 3+ ions statistically present inside the NCs 
will become trapped at the Si NC/Si02 boundary. These Er 3+ ions are susceptible of an 
instantaneous excitation directly into the first excited state, via direct absorption of photons 
with enough energy to excite erbium and create an electron-hole pair in the NC. In silicon 
the absorption processes are usually accompanied by phonon emission due to indirect band 
structure. When an Er 3+ ion is located inside a NC the absorption can accompanied by 
an excitation of this ion instead of phonon emission. This process is not considered in the 
theoretical part of this paper. We also note that the Er 3+ ions remaining in the NC will 
induce a donor center, as a result of which they will de-excite non-radiatively very fast j^] . 

As we can see from Fig. [H in panel b, the fast excitation of Er 3+ ions is completed 
within 20 ns after the excitation pulse, and we have estimated in our previous contribution 



2l| that ~ 50% of the Er 3+ content is involved in this fast process. The observed fast 
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(sub-microsecond) Er-related PL is only possible when Er 3+ is excited directly into the first 
excited state 4 Ii3/2- From Fig. [12] one can see that approximately 10 % of Er 3+ ions occur 
at a distance less than 1.15R from the NCs' centers, where according Fig. [TT] the most 
effective carrier cooling processes can provide Er 3+ excitation at times shorter than 20 ns. 
However, we should point out that the presented calculations of transition probabilities serve 
mostly as demonstration that such a fast mechanism is physically feasible. The theoretical 
considerations are valid for a single hole/single electron transition and for strictly spherical 
Si quantum dots. In this case, strong selection rules appear and only matrix elements 
of some higher multipoles of the electric field potential lead to considerable non-vanishing 
contributions, causing a strong decrease of the probability with increasing distance from Er 3+ 
to Si NCs. E.g., the Auger process accompanied by the 3e± — > 2eo transition is dominated by 
the quadrupole-dipole interaction (the dipole is related to Er 3+ ). Therefore one expects that 
the real probabilities might be higher than the calculated ones, and will decrease less abruptly 
with distance, due to non-sphericity of Si NCs. Even more important is the possibility of 
multiple exciton generation in a single NC, also not accounted for in the calculations. Under 
the excitation conditions used in the experiment, we estimate that on average about 5 
electron-hole pairs are created per Si NC. On one side it will directly increase the excitation 
probability of Er 3+ by a given NC at the initial stage. On the other side the strong exciton- 
exciton interaction will then generate very hot carriers leading to population of higher excited 
states whose participation has not been considered in Section IV. This will result in much 
higher transfer rates, due to a higher electron state localization outside Si NC. All these 
effects will contribute to the "fast" PL and improve the agreement with our experimental 
results. 

Previously, in 21|, we have estimated that ~ 0.5% of the total Er 3+ decay radiatively 
(r ~ 3 ms) after excitation. From these 0.5% of ions, ~ 50% have been excited due to 
the fast process, but have not been de-excited by the inverse mechanism: the "residual 
concentration". The other ~ 50% (~ 0.25% of the total Er 3+ content in the sample) have 
been excited due to inter-band recombination of the carriers in the NC. In our theoretical 
considerations we have shown that the excitation via the Auger process, accompanied by 
recombination of an exciton in the Si NC, can only take place via the contact. The Forster 
mechanism is not effective as shown in Sec. IVC. In any case the excitation of Er 3+ ion 
accompanied by a recombination of an electron-hole pair takes place into the second or 
the third excited state of the 4f shell of Er 3+ . The characteristic excitation of the 1.5 fim 
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Er 3+ PL is, in this case, a two-step process with time constant r = T\ + T2, where t\ is 
the characteristic Er excitation into one of the higher excited states by the band-to-band 
recombination of confined carriers, on the nanosecond or even shorter time scale, and Ti ~ 
1.2 [is is the relaxation time of the excited Er 3+ from a higher into the first excited state 
4 Ii3/2- This process is responsible for the fis rise of the Er-related PL signal. 

Following the described model, one expects that the fast-to-slow Er 3+ PL ratio will depend 
on the number of "hot" carriers confined in the Si NCs: a high number of "hot" carriers 
will favor the fast intra-band excitation process and thus increase the relative importance 
of the "fast" Er 3+ PL. This is indeed confirmed in Fig. El where we can see that under 
excitation conditions where creation of "hot" carriers is more likely, i.e. excitation with 
higher energy photons or high flux pumping and subsequent multiple carrier generation in 
the NCs, the fast-to-slow PL intensity ratio of Er PL increases. In the case of formation 
of several electron-hole pairs per Si NC, these can undergo a quick Auger recombination 
process, and in that way the excess pairs escape from participating in the "slow" excitation 
transfer process. 

Therefore, the following sequence of excitation and de-excitation processes can be pro- 
posed: 

1. After the laser pulse, an Auger-like process of fast excitation takes place, by intra- 
band relaxation of "hot" carriers transferring their excess energy directly into the first 
excited state of Er 3+ ion. Up to ~ 50% of all the Er 3+ content is in the effective range 
of this interaction, before the de-excitation processes start to decrease the number of 
excited ions. We note that this perce ntag e is dependent on the erbium concentration 



as indeed observed experimentally in 39] 



2. Fast de-excitation of Er 3+ takes place by a reverse process, transferring energy to 
carriers confined in the NCs. This process must be phonon assisted when the involved 
transitions do not match the energy conservation requirements; which will manifest in a 
slower PL decay and temperature dependence of the characteristic decay constant. The 
"hot" carriers in the NCs undergo an intra-band relaxation accompanied by phonon 
emission, which reduces the number of carriers available for the "quick" excitation 
process. The rate of relaxation processes should increase with temperature. 

A small percentage of the excited Er 3+ ions will overcome the fast de-excitation, giving 
rise to the "residual concentration" of Er 3+ which is excited in the nanoseconds time 
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window, and de-excites radiatively. This is revealed in the Er 3+ PL kinetics as the 
non-zero origin of the microsecond rise and subsequent radiative decay of Er 3+ PL. 

3. Once the confined carriers have cooled down to the bottom (top) of the conduction 
(valence) band, there is a probability of Er 3+ excitation by recombination, to the 
upper Er excited states, as described above. This probability exists also on shorter 
time scales, but the number of Er 3+ ions which can be accessed is low, and therefore 
this inter-band process is covered by the intra-band excitation process. 

Based on the presented experimental data and theoretical modelling, we conclude that 
three types of optically active Er 3+ ions coexist in the investigated material: 

• Type 1: Er 3+ ions that can only be excited resonantly, under direct excitation. These 
ions have predominantly radiative decay — like Er 3+ in SiC>2 — and are at large enough 
distance from the Si NCs to prevent their interaction. 

• Type 2: Er 3+ ions that are excited by energy transfer from inter-band recombination 
of excitons in Si NCs, at the microseconds time regime and which decay predominantly 
radiatively (r ~ 2-3 ms), independent of temperature. This type of Er 3+ ions are those 
accounted for in the usual estimations of optical activity measurements and constitute 
about 0.25% to 0.5% of the total Er 3+ content. 

• Type 3: Er 3+ ions excited via the intra-band transition of carriers in the NCs and with 
very strong non-radiative quenching, responsible for the sub-/is PL described above, 
whose properties mirror those of Er 3+ in crystalline- Si (~ 50% of Er 3+ ). A small 
percentage of these ions, will overcome the fast non-radiative de-excitation, and will 
be also accounted for in optical activity measurements. 

Finally, for the sake of completeness we note that in addition to these, there could also 
be Er dopants which are not optically active due to, e.g. precipitation. 

An independent confirmation of the proposed excitation model of Er 3+ ions by intra-band 
transitions of confined carriers in the Si NC is indeed given by Fig. [HJ There, for excitation 
energy higher than E^ pa 2.6 - 2.7 eV, a second excitation mechanism (due to intraband 
carrier cooling) is enabled, increasing the ratio between absorption and effective excitation 
cross section of Er PL. This threshold energy is sufficient to create a "hot" carrier that can 
excite Er 3+ directly into the first excited state by cooling into the bottom of the conduction 
band (or the top of the valence band) Q|. 
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VI. CONCLUSIONS 



With the results of this study, an important puzzle concerning the mechanism of excita- 
tion of Er 3+ by Si NCs has been solved. We have shown that the "missing" dopants that 
were apparently losing optical activity upon doping with Si NCs, and which did not con- 
tribute to PL, are actually very efficiently excited by the Si NCs, via an intra-band Auger 
transfer process, but undergo also a very effective excitation back-transfer process. The 
back-transferred carriers can again excite the Er 3+ ion or escape from being available for the 
intra-band excitation process due to a thermalization and Auger recombination processes 
inside a NC. In particular, we point out that the results of this study clearly show that in the 
first microsecond after the excitation laser pulse, a vast majority of Er dopants attain the ex- 
cited state; that implies that in that short time window the population inversion is reached, 
in that way fulfilling a necessary although in itself insufficient condition for realization of op- 
tical gain and laser. Future research will tell, whether the Si02:(Er, Si NCs) material can be 
engineered in such a way that the sensitization of Er emission is realized without the detri- 
mental effect of reduction of the concentration of Er 3+ ions, with the temperature-stable, 
predominantly radiative recombination. This study shows that achieving that will require 
careful and simultaneous optimization of Er 3+ and NC concentrations, NC size and size 
distribution, and, very importantly, very precise tuning of Er - NC distance, on a nanometer 
scale. 

On the other hand, we point out that fast recombination kinetics reported here leads to 
the effective recombination rate of Er 3+ in the environment of Si NCs exceeding by 2 orders 
of magnitude the fastest quenching rate (due to the Auger process involving free carriers) 
reported, to our knowledge, for Er in crystalline Si. This very efficient PL quenching of 
Er-related PL might be explored for GHz modulation of the 1.5 /zm emission from Er-doped 
structures. 
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APPENDIX: CALCULATIONS FOR THE EXCITATION DUE TO INTRA- 
BAND TRANSITION 



When dealing with the Coulomb interaction between an /-electron of an Er 3+ ion situated 
in SiC>2 and a carrier confined in Si NC, one should take into account the difference in 
dielectric constant values of Si (e\ = 12) and SiC>2 (£2 = 2). We note that the Auger process 
is determined by the high frequency dielectric constant, as the transition energy Aff is 
much larger than the lattice vibration energy 



The potential created by a point charge q at the distance a from the center of the sphere 
of radius R (R < a) with dielectric constant E\ in the media with dielectric constant e<i is 
obtained as a solution to the Poisson equation and is given by the equations: 

Q 



$i(r,a) 

e 2 a 



^ \a) le x + (Z + l)e 2 



(A.l) 



1=1 

inside the sphere (r < R), and 

- 5Fnr - ^fi (f i^kw; (*)'«• (A ' 2) 

outside the sphere (r > R), where $ is the angle between r and a: cos($) = (r, a)/ra. 
Note that the Coulomb potential in the case of interaction of two charges is given by 
I [$(r, a) + $(a, r)], but it is easy to show that in our case $(r, a) = $(a, r). 

Let us consider the potential $ in Eqs. (EJ) , (JTJ) . Introducing the coordinate r' related to 
the center of the ion (r' 77, where 77 is the size of the /-shell) and using the fact that 
Tf <C a, the potential can be expanded into a series by r' taking into account the linear term 
only: 

$(r, a + r') w $(r, a) + ^l^ r '. (A.3) 

oa 

In order to use formulae ( 1A.1I) and (1A.2I) the integration of matrix elements in Eq. (JTj) should 
be produced over r' for < r' < rj. Relatively high energy barriers at the boundary of 
a NC (3.2 and 4.3 eV for electrons and holes, respectively) allow only a small portion of 
the confined carriers charge density to penetrate outside. The charge density of confined 



carriers occurring in Si02 due to tunnelling is just a few percents [22J. Therefore the largest 
contribution is given by the Coulomb interaction induced by the carrier density inside the 
NC. To that end, it is enough to use potential $1 (see Eq. lA.lj) only. We can write for it 
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where 



and 



Si r 

j = — Ji + - j 2 , 

a r 



(2/ + l) £2 /r 



^ Zei + (/ + l)e 2 \a 



(I + l)P;(cos^) + cos$ 



0P,(cos#) 
<9cos"$ 



(2Z + l)e 2 /ry aFl(costf) 



Zei + (/ + l)e 2 ^a/ <9cos-# 



(A.5) 

(A.6) 
(A.7) 



Then Eq. (J7J) transforms into 



2ire 4 1 



a 1 1 M,M' ' ff 



E ^Ei d //'^ M 'i J i^)i 2 j r (iv), 



(A.8) 



where the ion dipole momentum is given by 



d ff = j ^(r)Wv(r)d 3 r. 



(A.9) 



Averaging Eq. flA.8[) over the directions of dffi one gets 



(A.10) 



where the dimensionless factor /^(a) is defined by 



= 4rf-) E Ki'.M'lJKM)] 5 



A",: V a 



(A.H) 



M.M' 



where the square of the matrix element absolute value is assumed to be averaged over the 
directions of vector a. Introducing the radiative lifetime r ra( j of the erbium ion in the first 
excited state ( 4 Ii3/ 2 ) [41 ] 

1 _ 1 2,2 V^( A //'' 3 



Trad Nf, 3 



ft 4 c 3 



(A.12) 



one gets finally Eq. (Q. The results of calculations of the factors ij'j(a) for electrons and 
holes are given in Fig. [101 

For completeness, we consider the contribution of carrier density outside the dot to the 
probability of Er 3+ excitation. The potential <J> 2 given by Eq. ( 1 A. 21) should be used outside 
the NC. We get 



<9$ 2 (r,a) q 



da 



R 2 e 2 



(A.13) 
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where 



_ aP , rP , 

" 9 1 J 2 

a 2 ar 



a)P 2 



r — a 



£1 - £2 



let + l)e 2 



Jo 



p\ 

a) 



(/ + l)P(costf) + costf- 



£1 — e 2 /T y <9p (cos 1?) 

dcos?? 



dcosfi 



(A.14) 
(A.15) 

(A.16) 



/Jl ^ + + 1 ) 5 2 W 

Producing the calculations analogous to the ones described above for the carriers being 
inside the NC, one can find that the contribution of the confined carriers tunnelling to the 
excitation probability is given by an expression similar to Eq. flA.lOj) . Our calculations have 
shown that the input is negligible. 



27 



[1] A. J. Kenyon, P. F. Trwoga, M. Federighi, and C. W. Pitt, J. Phys.: Condens. Matter 6, L319 
(1994) 

[2] M Fujii, M. Yoshida, Y. Kanzawa, S. Hayashi, and K. Yamamoto, Appl. Phys. Lett. 71, 1198 
(1997) 

[3] C. E. Chryssou, A. J. Kenyon, T. S. Iwayama, C. W. Pitt, and D. E. Hole, Appl. Phys. Lett. 
75 2011 (1999) 

[4] D. Pacifici, G. Franzo, F. Priolo, F. Iacona, and L. Dal Negro Phys. Rev. B 67, 245301 (2003) 

[5] P. G. Kik, M. L. Brongersma, and A. Polman, Appl. Phys. Lett. 76, 2325 (2000) 

[6] D. Kovalev, H. Heckler, M. Ben-Chorin, G. Polisski, M. Schwartzkopff, and F. Koch, Phys. 

Rev. Lett. 81, 2803 (1995) 
[7] M. E. Castagna, S. Coffa, M. Monaco, L. Caristia, A. Messina, R. Mangano, and C. Bongiorno, 

Physica E 16, 547 (2003) 
[8] A. Nazarov, J. M. Sun, W. Skorupa, R. A. Yankov, I. N. Osiyuk, I. P. Tjagulskii, V. S. 

Lysenko, and T. Gebel, Appl. Phys. Lett. 86, 151914 (2005) 
[9] R. J. Walters, G. I. Bourianoff, and H. A. Atwater, Nature Mat. 4, 143 (2005) 
[10] T. Forster, Ann. Physik 2, 55 (1948) 

[11] M. J. A. de Dood, J. Knoester, A. Tip, and A. Polman, Phys. Rev. B 71, 115102 (2005) 
[12] V. M. Agranovich and M. D. Galanin, Electronic Excitation Energy Transfer in Condensed 
Matter, Elsevier (1982) 

[13] X. L. Wu, Y. F. Mei, G. G. Siu, K. L. Wong, K. Moulding, M. J. Stokes, C. L. Fu, and X. M. 

Bao, Phys. Rev. Lett. 86, 3000 (2001) 
[14] R. A. Senter, C. Pantea, Y. Wang, H. Liu, T. W. Zerda, and J. L. Coffer, Phys. Rev. Lett. 

93, 175502 (2004) 

[15] M. Fujii, K. Imakita, K. Watanabe, and S. Hayashi, J. Apl. Phys. 95, 272 (2004) 

[16] O. Savchyn, F. R. Ruhge, P. Kik, R. M. Todi, K. R. Coffey, H. Nukala, and H. Heinrich, Phys. 

Rev. B 76, 195419 (2007) 
[17] K. Imakita, M. Fujii, and S. Hayashi, Phys. Rev. B 71, 193301 (2005) 
[18] J. Lee, J. H. Shin, and N. Park, J. Lightwave Technol. 23, 19 (2005) 

[19] M. Wojdak, M. Klik, M. Forcales, O. B. Gusev, T. Gregorkiewicz, D. Pacifici, G. Franzo, F. 
Priolo, and F. Iacona, Phys. Rev. B 69, 233315 (2004) 

28 



[20] P. G. Kik and A. Polman, J. Appl. Phys. 88, 1992 (2000) 

[21] I. Izeddin, A. S. Moskalenko, I. N. Yassievich, M. Fujii, and T. Gregorkiewicz, Phys. Rev. 
Lett. 97, 207401 (2006) 

[22] A. S. Moskalenko, J. Berakdar, A. A. Prokofiev, and I. N. Yassievich, Phys. Rev. B 76, 085427 
(2007) 

[23] A. J. Kenyon, C. E. Chryssou, C. W. Pitt, T. Shimizu-Iwayama, D. E. Hole, N. Sharma, and 

J. Humphreys, J. Appl. Phys. 91, 367 (2002) 
[24] A. Schenk and G. Heiser, J. Appl. Phys. 81, 7900 (1997) 
[25] J. R. Chelikowsky and M. Schlueter, Phys. Rev. B 15, 4020 (1977) 

[26] V. N. Abakumov, V. I. Perel, and I. N. Yassievich, Nonradiative Recombination in Semicon- 
ductors, Elsevier (1991) 

[27] E. Hendry, M. Koeberg, F. Wang, H. Zhang, C. de Mello Donega, D. Vanmaekelbergh, and 

M. Bonn, Phys. Rev. Lett. 96, 057408 (2006) 
[28] B. K. Ridley, Quantum Processes in Semiconductors, (Clarendon Press, Oxford) (1982) 
[29] A. A. Prokofiev, A. S. Moskalenko, and I. N. Yassievich, Mater. Sci. and Eng. B 146, 121 

(2008) 

[30] M. D. Shinn, W. A. Sibley, M. G. Drexhage, and R. N. Brown, Phys. Rev. B 27, 6635 (1983) 
[31] A. A. Prokofiev et al., J. of Lum. 121, 222 (2006) 

[32] A. A. Prokofiev, I. N. Yassievich, H. Vrielinck, and T. Gregorkiewicz, Phys. Rev. B 72, 045214 
(2005) 

[33] I. N. Yassievich, and L. C. Kimerling, Semicond. Sci. Technol. 8, 718 (1993) 

[34] A. S. Moskalenko, I. N. Yassievich, M. Forcales, M.Klik, and T. Gregorkiewicz, Phys. Rev. B 

70, 155201 (2004) 
[35] M. Cardona and F.H. Pollak, Phys. Rev. 142, 530 (1966) 

[36] A. Thraendhardt, C. Ell, G. Khitrova, and H. M. Gibbs, Phys. Rev. B 65, 035327 (2002) 

[37] B. R. Judd, Phys. Rev. 127, 750 (1962) 

[38] W. J. Miniscalco, J. Lightwave Techn. 9, 234 (1991) 

[39] K. Imakita, M. Fujii, and S. Hayashi, Eur. Phys. J. D 34, 161 (2005) 

[40] D. Timmerman, I. Izeddin, P. Stallinga, I.N. Yassievich, and T. Gregorkiewicz, Nature Pho- 
tonics 2, 105 (2008). 
[41] A. Suchocki, and J. M. Langer, Phys. Rev. B 39, 7905 (1989) 



29 



30 



FIGURE CAPTIONS 



FIG. 1: Time-integrated PL spectra at T — 10 and 300 K, and excitation wavelength X exc — 450 
nm. Si NCs, Er-related, and a third (at low T) PL bands can be observed. The Si NC-related 
band is blown up for clarity. The arrows show PL peaks related to emission from higher excited 
states of the Er 3+ ion, superimposed to the NCs excitonic PL. 

FIG. 2: Time-resolved contour plot of the PL spectrum of the Si NCs excitonic-related band, in 
the microseconds (a) and in the sub-microsecond (b) time ranges. The PL intensity (color contour 
plots) is represented as a function of time (X axis) and detection wavelength (Y axis). Panel a also 
shows the spectrum recorded at t = 100 /is, and the PL decay at A = 850 nm (sections from the 
contour plot). Drift of the band maximum toward longer wavelength is observed for longer time 
delay in panel a. Note the difference between the center of the bands in both figures (~ 200 meV) . 

FIG. 3: Room temperature PL decay kinetics recorded at A = 1200 nm (maximum of the PL broad 
band shown in Fig. [I]) and A = 1535 nm. In the inset, the spectrum recorded at t = 200 ns after 
the laser excitation pulse, at T = 10 K, is shown. 

FIG. 4: In panel a: Si NCs PL kinetics, recorded at A = 860 nm, and room temperature; time 
resolution of the system was 2 ns. Excitation was provided with a 5 ns pulse, at \ exc = 450 nm. 
Initial nanoseconds decay is stretched and followed by a final decay with characteristic lifetime of 
thick microseconds. In panel b: Er-related 1.5 /im PL kinetics for the first microsecond after the 
excitation pulse is shown, at T = 10 K and RT. 



31 



FIG. 5: Er-related PL kinetics, in double logarithmic scale for different excitation conditions, at 
RT. Intensity has been normalized to the maximum of the 'slow' component, in order to compare 
the fast-to-slow intensity ratio. In panels (a) and (b) the excitation wavelength (450 nm) has been 
kept constant, and flux has been decreased about one order of magnitude. In panels (b) and (c), 
flux has been kept in the same order of magnitude, and wavelength excitation has been increased 
to 650 nm. 

FIG. 6: Flux dependence of Er 3+ and Si NCs PL for several excitation wavelengths is shown, 
recorded at RT. 

FIG. 7: PL excitation cross section of Er 3+ and Si NCs PL as a function of excitation wavelength, 
which have been determined by using formula dU) to fit the curves of Fig. [6j 

FIG. 8: Er- and Si NC-related effective excitation cross section api as a function of the optical 
absorption coefficient a, for each given excitation wavelength. 

FIG. 9: Electrons and holes energy levels in Si NCs in Si02, as a function of NC diameter. On 
the right-hand side the Er 3+ energy levels are presented. The most effective erbium excitation and 
de-excitation processes due to intra-band carrier transitions are shown. 

FIG. 10: Factors 1^ for erbium excitation/de-excitation by electron intra-band transitions as a 
function of a/R, where a is the distance of the Er 3+ ion from the center of the NC and R the radius 
of the NC. 

FIG. 11: The probability of the most effective excitation and de-excitation processes as a function 
of a/R (a is the erbium distance from the center of a NC and R the radius of the NC) for a NC 
with diameter 3.1 nm. 

FIG. 12: A simulation of the distance distribution of Er 3+ ions to the center of their nearest 
neighboring NC assuming a random distribution of both: 6.6% of the Er 3+ ions are contained 
inside NCs. 
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